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Abstract We present a new dispersive framework for the extraction of the strong 
coupling constant a s from r-lepton decays. A new feature of our procedure is the 
use of the quark- hadron duality on the limited region s<j < s < m 2 . The dual- 
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ity point s<j and the MS strong coupling constant a s (m 2 ) are self-consistently ex- 
tracted from the r data for the non-strange vector spectral function. We use 2005 
ALEPH and 1998 OPAL experimental data on the vector spectral function. We com- 
pare the new framework with the contour improved perturbation theory up to or- 
der Qg. The new procedure yields systematically lower values for a s . From the 2005 
ALEPH data, we obtain a s (m 2 ) = 0.308 ± 0.014 eX p ± 0.005th which corresponds to 
a s (Af 2 ) = 0.1170 ± 0.0018 oxp ± 0.0007 th ± 0.0005 ov . The extracted value for the 
duality point is found surprisingly stable against perturbation theory corrections 
s d = 1.71 ± 0.05cx P ± 0.00 th GeV 2 . From the 1998 OPAL data, we obtain « s (m 2 ) = 
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0.290 ± 0.023 cxp and s d = 1.68 ± 0.10 cxp GeV 2 . 



1 Introduction 



■ Keywords tau lepton decay • renormalization group equation • perturbation theory 

" data analysis 

in 
Ov 
in 

(ST As is well known, the inclusive hadronic decays of the r-lepton may be reliable studied 

within perturbative QCD (see seminal work [I] and the literature therein). A general 
approach to analyzing the perturbative and non-perturbative aspects of the r-system 
observables is the renormalization group improved perturbation theory augmented with 
the Wilson's Operator Product Expansion (OPE) [2J. The characteristic energy scale 
is relatively small, m T ~ 1.8 GeV (m T being the mass of the r-lepton). Hence, the non- 
perturbative effects of QCD cannot be completely ignored. The original study [TJ has 
shown that they are small and under control within the OPE. In the following years, 
the inclusive hadronic quantities of the r system have been intensively exploited to 
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precisely determine the strong coupling constant a s (m T ). This became feasible because 
the observables of the r system are sensitive to the concrete value of a s and the accuracy 
of the experimental data for a variety of the observables has been considerably improved 
(for recent review see [3])- 

In past few years, the determination of the strong coupling constant from non- 
strange hadronic r-data has received a renewed interest. It was pointed out [4] that 
there is not good agreement between recent two highest precision low-energy deter- 
minations of a s . These determinations come from the finite energy sum rule (FESR) 
analysis of hadronic r decay data [5] and from a lattice perturbation theory analysis 
of ultraviolet-sensitive lattice observables [6]. 

a s {Ml) = 0.1212 ± 0.0011 (r decay) (1) 
a s (Ml) = 0.1170 ±0.0012 (lattice). (2) 

Moreover, different determinations of a s from the same r data [73[8][9] are not fully 
consistent within their mutual errors (see work [ID] and the literature therein). This 
discrepancy has stimulated a number of new theoretical investigations on the appli- 
cation of the FESR in r decays (see works [4lf5ll 10U 1 1U12II13U14II15U16U 17U 18U 19j . The 
standard FESR technique based on the truncated OPE series has been reconsidered. 
The small but still significant non-perturbative effects have been included into analysis 
[19] - On the one hand, the impact of the higher order terms of the OPE (neglected 
in the standard analyzes) has been estimated [4, 15, 16, 18 i. It was confirmed that their 
influence on the extracted value of a s is not small in the separate vector and axial 
vector channels. To suppress these contributions in the FESR the so-called pinched 
weights introduced 1.15 Hi lis]. On the second hand, using the physically motivated 
model [1011111119] the impact of the non-perturbative corrections coming from the pos- 
sible duality violations (DVs) [20] has been estimated. In the separate vector and axial 
vector channels the DVs was found to be appreciable (see recent work [15] and refer- 
ences therein). The pinched weights have also been employed to reduce the effects of 
DVs 18 . Possible non-perturbative corrections to the FESR (direct instantons, duality 
violation and tachyonic gluon mass) which cannot be described within the OPE have 
been estimated in |13| . 

As is well known, in the time-like region the renormalization group (RG) invariance 
cannot be used unambiguously. For this reason, two different methods are used to 
perform the RG resummation within the FESR. These are fixed order perturbation 
theory (FOPT) and contour improved perturbation theory (CIPT) |21II22| . These two 
approaches lead to differing results. The values of a s extracted from r decays employing 
CIPT have always been higher. A critical comparison of these two approaches may be 
found in recent works [231124] . In [M] FOPT was approved as a better approximation 
to the true result. In contrast, authors of [3] and [5] favored CIPT. 

Note that the non-physical singularities of the perturbative running coupling (the 
Landau pole problem) which occur at small space-like momenta may, supposedly, de- 
teriorate the extracted values of the parameters [25]. In particular, CIPT suffers from 
this shortcoming |24| . To cope with this problem dispersive or analytic approaches to 
perturbative QCD have been developed |2^|2Tl[2^pgpJl|3TT|3^^ . 
In works [27] and [28] , the r lepton decay rate has been analyzed within a simple and ef- 
fective dispersive technique, Analytic Perturbation Theory (APT) (for reviews see |29l 
38,39 ). However, the minimal analytic QCD model (the same APT) predicts, from 
the non-strange r lepton decay data, too large value for the strong coupling constant, 
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a s (m^) = 0.403±0.015 [28] . The advantages and shortcomings of the three approaches 
to the r decays (FOPT, CIPT and APT) were thoroughly analyzed in 730] • It should 
be noted that APT as well as its generalized versions suggested more later 31,33,34, 
35,36,37 proved to be very useful from the phenomenological point of view. A re- 
markable feature of these modified expansions is the better convergence and improved 
stability property with respect to change of the renormalization scheme. Nevertheless, 
one should keep in mind that an analytic approach based only on perturbation theory 
can not be defined unambiguously, in fact, there is not a unique recipe for removing 
the Landau singularities from the running coupling. 

In our earlier work [40] we have suggested a dispersive approach to analyze the r 
decay data. In contrast to CIPT, the new approach is based on the improved approxi- 
mations to the Adler function which incorporate correct analyticity and RG invariance 
properties of the exact function. Moreover, the approximations correctly reproduce 
the required ultraviolet and infrared properties of the exact Adler function. Another 
feature of the new framework is the use of the quark-hadron duality in the limited 
region < s < m T . The QCD scale parameter ^%jg and the duality point s d may be 
determined, self-consistently, from the experimental data |40| . 

In the present article, we investigate the new framework more thoroughly. We revise 
part of the results of work [40] . We present a more accurate test of the convergence of 
the numerical results in perturbation theory. The numerical value of the duality point 
is found to be remarkable stable with respect to higher order QCD corrections. More 
importantly, we study the stability of the results with respect to small change of the 
experimental data. In Sect. 2 we critically analyze the FOPT and CIPT approaches 
to the r-decay. A dispersive modification of the CIPT suggested in [3D] is discussed in 
more detail. In Sect. 3 we give corrected numerical values for a s and s<j extracted from 
the 2005 ALEPH data. We thoroughly investigate the stability of the results comparing 
the new and CIPT determinations of a s order by order in perturbation theory. In addi- 
tion, we analyze the ALEPH non-strange data employing the renormalization scheme 
invariant (RSI) framework suggested in [41] . We also analyze 1998 OPAL [9] vector 
data within the new dispersive framework. Conclusions are summarized in Sect. 4. 



2 Theoretical Framework 

Let us briefly recall some basic facts about the QCD analysis of the hadronic decays 
of the r-lepton through the FESR [42]. The non-strange vector component of the r- 
hadronic width is determined as 

2 f m * 

Rr,v = 6\V ud \ Sew / w T (s)vi(s)ds, (3) 
Jo 

where 

V u( \ is the flavor CKM matrix element, Sew denotes a short-distance electroweak 
correction Q and v\{s) is the vector spectral function defined through the correlation 



In what follows, we neglect the small additive electroweak correction 8' ew- 



1 



function 

vi(s) = 2nlmn udy (-s). (4) 

It is more convenient to define a renormalization scale invariant quantity, the Adler 
function 

D(Q 2 )^-47v 2 Q 2 ^n ad , v (Q 2 ), (5) 

here, we have defined s = q 2 = — Q 2 . In the exact theory, the correlation function 
n a( i y(z) and the Adler function are analytic functions except the cut running along 
negative 2-axis. This implies the FESR relation 

^ y =-£i^ EW £^ e (i-^) (i+ (6) 

here, the integration contour is a circle of radius so (so = m 2 ). In the case of massless 
quarks, the Adler function has the perturbation theory expansion [23] 

2 °° 2 " +1 1 Q 2 

D{Q 2 ) = a s (pL 2 ) n V* kc nk L _1 where L = In — , (7) 
™=0 fc=l 

a s (n 2 ) = a s (/i 2 )/7r and a s (A* 2 ) denotes the strong coupling constant normalized at the 
scale /j,. It follows from the renormalization scale invariance of the Adler function that 
only the coefficients c„i are independent. All other coefficients are determined in terms 
of the c„i and /3-function coefficients through the RG equation [231124] . In practice the 
series (0 is truncated at some finite order. 

The approximations to the Adler function obtained by truncation of the series (0 
have correct analytical properties of the exact function. In the case of FOPT, the series 
(0) is inserted into contour integral (01 and integrated term-by-term. Afterwards, the 
normalization scale is determined choosing /i = m T [23] ■ However, we could start from 
the original formula ([3]) with perturbation theory expansion for the spectral function. 
The expansion for the spectral function is obtained by insertion series (0) into inversion 
formula ()15|) (see below) and integrating term-by-term. So, we could achieve the same 
result without using the FESR relation (0. Thus, within FOPT, formulas (0) and (01 
are equivalent. However, the approximations to the Adler function employed within 
FOPT do not describe correctly the asymptotic behavior of the exact function for 
Q 2 —> oo. In the standard analysis of the r data this fact is irrelevant. However, as we 
shall see latter, this is not the case for the new framework accepted in this paper. 

Since the Adler function is the renormalization scale invariant quantity, one may 
choose /i 2 = Q 2 in series (0. Thus, one obtains the RG improved expansion 

oo 

D(Q 2 )\bg = 1 + d(Q 2 )\ BG = 1 + £ d ^(Q 2 ) (8) 

n=l 

where d n — c„i, a s (Q 2 ) = a s (Q 2 )/n, a s (Q 2 ) being the running coupling. The first two 
coefficients in the expansion (0) are the renormalization scheme invariant. The known 
coefficients in the MS scheme for nj = 3 quark flavors take values d\ = 1, di ~ 1.6398, 
c?3 ~ 6.3710 and ~ 49.0757. The last coefficient was calculated recently by the 

2 We use the normalization of the spectral function with the naive parton prediction 

^1, naive ~ 1/2. 
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authors of |12] by using powerful computational techniques. The approximations to 
the Adler function constructed by truncation the series j8} have correct ultraviolet 
asymptotical behavior (d(Q 2 ) — > as Q 2 — > oo), however they violate the cut-plane 
analyticity of the exact Adler function due to the non-physical Landau singularities of 
the perturbative running coupling. One may assume, without loss of generality, that the 
running coupling has only one Landau singularity located on the positive Q 2 axis [40] . 
This is true for the asymptotic solutions and for more accurate Lambert-W solutions 
to the RG equation [431144 ||4"5] . Then we may derive (see [40] ) the violated dispersion 
relation for the QCD correction to the Adler function: 

d(Q 2 )\ RG = d(Q 2 )\ APT +d L (Q 2 ), (9) 

where the function d(Q 2 )|APT satisfies the normal DR 

rf(Q 2 )lAPT = - r ££ ^da, (10) 

t Jo ° + Q 

with the effective spectral density determined as 

PcffM =Im{d(-cr-je)| R G}- (H) 
It is to be noted here that the function 

J D(Q 2 )APT = l + rf(0 2 )lAPT (12) 

is the analytic image of the perturbative Adler function determined in the sense of 
the Analytic Perturbation Theory (APT) approach of Shirkov and Solovtsov [261127] . 
The second term in (J9]) is the contribution coming from the Landau singularity. It is 
represented by the contour integral [40] 

here, the integral is taken round the circle {C : C — S L + S L ex P WO) ~ n < <j) < ir} in 
the positive (anti-clockwise) direction, with sl being the Landau singular point. 

In the popular framework, referred to as contour improved perturbation theory 
(CIPT) [21122] , the (truncated) expansion © is inserted into the FESR integral © 
and then integrated term by term. At this point, one ignores the fact that with the 
approximation ([S]) formulas © and are not equivalent. Indeed, the FESR relation 
© can not be derived because of violated analytical properties of the approximation 
([8]). This inadequacy repeatedly discussed in the literature (see for example [24] and 
30 ). Nevertheless, CIPT has been very successful from the phenomenological point 
of view. On the other hand, APT is free from this drawback. However, the analysis 
of the r decay data based on APT with massless quarks gave too large value for the 
strong coupling constant |28| . Furthermore, in the infrared region, the Adler function 
can not be reproduced correctly within APT, CIPT or FOPT. Thus, in APT the 
running coupling a s (Q 2 ) has a finite limit as Q 2 — > [26] . This leads to the apparent 
contradiction in the case of the Adler function. In fact, the Adler function should vanish 
at Q 2 = 0, as is manifested by Chiral Perturbation Theory [47]. In work [31], APT has 
been modified by considering the quark mass threshold effects for the light quarks. In 
this way, correct descriptions of the Adler function and r data was achieved. However, 
too large values for the effective quark masses (m u ~ ~ 330 MeV) was predicted. 
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Fig. 1 The integration contour in the complex Q 2 plane used in the inversion formula 115t 
in the case of the approximation £ ) (Q 2 )|rg violating the DR. Branch points on the real axis 
are represented by the blobs and branch cuts by the zigzagging lines, denotes the Landau 
singularity. 



As is well known, in the exact theory the Adler function satisfies the dispersion 
relation (DR) 



D(Q 2 ) = Q 2 



2vi(s)ds 

Jo (s + Q 2 ) 2 



(14) 



the corresponding inversion formula reads 



— S+2€ 



D(z) 



dz, 



(15) 



where the path of integration, connecting the points — s =p ze on the complex z-plane, 
avoids the cut running along the real negative z axis. The integral being traversed in 
a positive (anticlockwise) sense. From the violated DR we may also derive the 
integral representation 



R.G 



s)ds 



(s + Q 2 ) 2 



(16) 



where the singular integral at the lower bound should be treated in the sense of dis- 
tribution theory 0. It is to be noted that the spectral function Vi^ G (s) may be again 
calculated via the inversion formula (|15|l , but now the integration contour should also 
avoid the non-physical cut running along the positive interval < z < (see Fig. 1). 

The dispersion relation ()14[1 may be used to construct the approximations to the 
Adler function with correct analyticity properties. To approximate the hadronic spec- 
tral function, one may use the global duality ansatz employed previously in works 46, 

H3 



«i(a) = 6(s d - s)v^(s) + 9(8 - Sd K QCD (s) 



(17) 



where w^ CD (s) is the perturbation theory approximation to the spectral function, 
«" p (s) denotes the non-perturbative component of the spectral function confined, pre- 
sumably, in the low energy region, and is the onset of perturbative continuum, an 



3 We have confirmed formula JT6j in the case of the one-loop order /3-function. 
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infrared boundary in Minkowski region above which we trust pQCD 0. One may also 
construct a "semi-experimental" spectral function 

v^"(s) = 9(s d - s)v^(s) + 6(8 - Sd K QCD ( S ), (18) 

where v^ xp (s) denotes the genuine experimental part of the total "semi-experimental" 
spectral function. It was measured with high precision by ALEPH 08] and OPAL 
[9] collaborations in the range < \fs < m T = 1.777 GeV. Formula (|18[) extends 
the spectral function beyond the range accessible in the experiment. Formulas (|17|1 
and (I18|) provide practical realizations of the concept of the quark-hadron duality 
(see the original work [36]). In [47], this ansatz was used to determine the duality 
point s d for a given value of ^%jg, the QCD scale parameter in the MS scheme. The 
perturbative component ii^ CD (s) was constructed from the FOPT series (O choosing 
the normalization scale fi 2 = s d . Such a framework may be considered as a modification 
of FOPT. In [3D], we have used the same ansatz for the spectral function. However, our 
strategy was somewhat different. Starting from the ansatz (|18[1 . we have determined 
the parameters A^jg- and s d self-consistently from the r data. In contrast to [47] . 
we have used the RG improved approximation vf" G (s) to the spectral function. The 
function ff" G (s) is calculated by insertion of the (truncated) RG improved series ^ 
into inversion formula (|15l) . For s > 0, one finds [40] 

pQCD, s RG / \ APT/ n /, n \ 

«i («)=«1 0) = "l (s), (19) 

where d^ PT (s) is the spectral function determined in the sense of the Shirkov-Solovtsov 
APT x 

v£ PT (s) = -{l + r(s)) where r(s) = -[ PcS ^ da. (20) 
2 T J a G 

It follows from the duality relation (I18|l that one may calculate in QCD perturbation 
theory the decay rate of the r lepton into hadrons of invariant mass larger than y/s^ 



so that 



R^ D \s>s d =(>\V ud \ 2 S EW w T (s)vt PT (s)ds,^R c T J\ s>Sd (21) 

J s d 

$t{s a ,A 2 )= w T {s)vt PT (s)ds= w T (s)v c 1 xp (s)ds. (22) 

J s d J s d 

Using relation (|20p . one may express the left hand side of (1221) in terms of the effective 
spectral density [3D] 

# T (, d ,^)^(i^g d) 3 (i+g d ) (i+ : (sd)) 



47T 



1 -1, 3 2 (23) 

y (l-J/) (1 + v)Pce{m T y)dy, 



where s d = s d /m 2 . 

Inserting the duality ansatz (I18p into DR (I14|l one constructs the "semi-experimental" 
Adler function 

£>(Q 2 )Kcx P " = J D(Q 2 ,s d )| CX p + D(Q 2 ,s d )lpQCD, (24) 



4 It is assumed that < s c \ < m^. 
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where the experimental and QCD components of the Adler function are determined by 



D(Q 2 ,s d )\ cxp 



, [ a * 2v c 1 xp (s)ds ,2m 2 f°° 2wf QCD (s)ds 

Jo W ' Sd) ^ CD=0 ./ ( S + Q2)2 ' 



(25) 

In general, the QCD component uP may contain the non-perturbative correc- 

tions coming from the OPE as well as the duality violating terms [19] not included into 
the OPE. Intuitively, it seems to us that the non-perturbative corrections are more es- 
sential in the region < s < s d . In what follows, we will ignore these non-perturbative 
corrections into QCD component of the spectral function and employ the perturba- 
tive approximation (jT9j) . The power suppressed part of the "semi-experimental" Adler 
function is defined as 

D(Q 2 )\ pw . s = D(Q 2 )\< %exp „ - D(Q 2 )\ RG , (26) 
it may be represented in the form 40 

£(Q 2 )Ipw. s = 2 K(Q 2 ,s)(v c * p (s) - vt PT (s))ds - d L (Q 2 ), (27) 
Jo 



where K(Q 2 ,s) = Q 2 /(Q 2 + s) 2 . Formula (|27l) enables us to derive the asymptot 



it- 



expansion at large Q 2 



D(Q 2 )\ pw . s ~ Yl V* ( tS ) ( 28 ) 



where A = A^m is the QCD scale parameter in the MS scheme and the coefficients 7/ n 
depend on the dimensionless ratios A 2 /m 2 and s d /m 2 . In the case of massless quarks, 
the gauge invariant operator of dimension two can not be constructed. Hence, it follows 
that rji = 0. This condition from the OPE leads to the equation relating the parameters 
S(j and A with the experimental spectral function [3D] 

<P as (s d ,A 2 ) = ± T v c * p (s)ds, (29) 
m4 Jo 

where 

<Ms d ,^ 2 ) = %l + r( Sd )) + -! T r Posada +^^, (30) 

2 2-KITlr Jq 2 rtlr 

with s d = Sd/m 2 and the coefficient cl is a positive number independent of A 

c L =A- 2 ^-(f d(C)| RG rfC=^-% / d(sL + s L e^)| RG d^ (31) 

here sl being the Landau singularity of the running coupling. It is proportional to 
A 2 E 

Numerical values of the coefficient cl calculated in the MS scheme are listed in 
Table [1] In the calculations we have used the approximations to the Adler function of 
increasing order 0. All approximations have been constructed with the four- loop order 



5 Analytic expressions for sl in the MS scheme up to fourth order in perturbation theory 
may be found in 1451 . 

6 We use the abbreviation N^LO to denote the order O(oj^ 1 ) approximation to the Adler 
function. 
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Table 1 Numerical values of the coefficient cl in the MS scheme calculated with the four-loop 
order exact numeric running coupling. 

Approximations to the Adler function 

LO NLO N 2 LO N 3 LO N 4 LO 

c L 0.301262 0.453421 0.555401 0.651373 0.721687 



exact (numeric) running coupling. For the unknown 0(a^) correction to the Adler 
function, we use the geometric estimate c?5 = d^di/ds) = 378 [5]- 

An important remark is in order here. The advantage of the approximation «^ PT (s) 
is that it correctly describes asymptotic behavior of the exact function as s — > oo; in this 
limit v 1 (s) —¥ 1/2. In contrast, the FOPT approximation vf (s) increases with 
s as a polynomial of In s. This shortcoming of FOPT is irrelevant as far as the duality 
relation (|22[) is concerned. However, Eq. (129[) depends on the ultraviolet properties of 
the Adler function. This discussion suggests that a more consistent framework should 
be constructed in the contour improved scheme. In this work, we will refer the new 
framework as dispersive contour improved perturbation theory (DCIPT) Q. Although 
technically DCIPT resembles APT, there are significant differences between the two 
frameworks. Thus, in DCIPT we do not mention modifications of the QCD /3-function 
and running coupling. 

3 Numerical Results 

The parameters s d and A may be extracted from the data by solving the system of 
equations 

$ T (s d ,A 2 )= / w T (s)v c 1 xp (s)ds, (32) 

Js d 

$ as (s d , A 2 ) = -L f " v c 1 xp (s)ds, (33) 

where the functions & r and $ a s are defined in formulas (|23|) and (|30p . The right hand 
sides of Eqs. (|32p ~ (|33p are determined in terms of the empirical function «° xp (s). We 
have reconstructed the experimental vector spectral function from the 2005 ALEPH 
spectral data for the vector invariant mass squared distribution 8 . This was done, with 
the values \V ud \ = 0.9746 ± 0.0006 and S E w = 1-0198 quoted in 0. To interpolate 
the spectral function between the fixed experimental values of the energy squared, we 
use cubic splines. Evidently, the mean values of the parameters should be determined 
from the mean value of «° xp (s). The error analysis is based on the system of equations 
(I32H - II33I ) 40] , To determine the experimental uncertainties on the extracted values of 
the parameters, we use covariance matrices provided by ALEPH. Unfortunately, in the 
earlier work [40], we used (inconsistently) the N 2 LO value cl = 0.555401 (see Table [TJ 
in all other orders. In this work, we present corrected results. 

In general, the system (|32p ~ (|33p has more than one solution. For phenomenological 
reasons, we look for the solution in the limited region 280 MeV < A < 420 MeV. In this 
region, the system has only one solution. In Table [2j we give the central values for the 



7 In | 40| . we used the abbreviation APT" 1 ". 
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Table 2 Central values for the parameters in the MS scheme extracted from the 2005 ALEPH 
vector t data order-by-order within DCIPT. These results correspond to the four-loop order 
running coupling. 

Observable Approximation to the Adler function 

LO NLO N 2 LO N 3 LO N 4 LO 

s d GeV 2 1.707 1.710 1.709 1.707 1.705 
A GeV 0.486 0.378 0.348 0.332 0.323 

a s {ml) 0.401 0.337 0.321 0.313 0.308 



Table 3 The changes of the leading term induced by the consecutive corrections in the series 
@U and J35]>. 



Perturbative orders 


NLO 


N 2 LO 


N a LO 


N 4 LO 


DCIPT 


15.9% 


4.0% 


2.2% 


1.2% 


CIPT 


19.6% 


4.7% 


2.7% 


1.4% 



parameters extracted from 2005 ALEPH data within the new (DCIPT) framework. 
Formally, we may write a series for the numerical value of the coupling constant as 
follows 

4 

Q s (m^)| N 4 LO = «s("V)Ilo + Ak > 

where = Q s (m^)| N kLo — a s(TO T )|Nk-iLo- Using the numbers listed in Table[2] we 
obtain the series 

a s {ml) |n?lo T = 0.401 - 0.064 - 0.016 - 0.009 - 0.005. (34) 
In [40], from the same data, we have obtained the CIPT series 

M™t)In4lo = 0- 485 ~ °- 095 ~ 0- 023 ^ °- 013 - 0.007. (35) 

In Table [31 we give the changes (in percents) of the leading term induced by the 
consecutive corrections in the DCIPT and CIPT series. One sees that the DCIPT 
series (|34[) converges more rapidly. 

In this paper, we will estimate only so called indicative theoretical errors. These 
are defined as a half of the last retained term in the series [41] . As pointed out in [41j . 
this definition of the error is heuristic and indicative. From the DCIPT series (I34|l . we 
obtain the estimates 

<Mnv)lNLO = 0.337 ± 0.016 cxp ±0.032 th 
a s (m?)j N 2 LO = 0.321 ± 0.016 cxp ± 0.008 th 
a s (m?)j N 3 LO = 0.313 ± 0.014 cxp ±0.004 th 

a s (m 2 T )\ N 4 LO = 0.308 ± 0.014 cxp ± 0.002 th ± (0.0045 d5 ), (36) 

here we have also included the experimental errors 0. In previous paper |40| . we have 
found from the same data in the case of CIPT 

«s(wv)!nlo = 0.390 ± O.Olloxp ± 0.048 th 
8 We use formulas for the error analysis derived in 1401 . 
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Table 4 Comparison of the DCIPT and CIPT r decay determinations of the strong coupling 
constant at the scale M% = 91.187 GeV. Two errors are given, the experimental (first number) 
and the error from the evolution procedure (second number). 



Approximation 


a s (Mf)| D ciPT 


a a (M z 2 )| CIPT 


N 2 LO 


0.1187 ±0.0019 ±0.0005 


0.1238 ± 0.0009 ± 0.0005 


N 3 LO 


0.1176 ±0.0018 ±0.0005 


0.1224 ±0.0009 ±0.0005 


N 4 LO 


0.1170 ±0.0018 ±0.0005 


0.1217 ±0.0009 ±0.0005 



a s (>4)| N 2Lo = 0.367 ± 0.009 oxp ± 0.012 th 
« S ("1t)|n3 L o = 0.354 ± 0.008ox P ± 0.007 th 

a s (m 2 T )\ N 4 LO = 0.347 ± 0.008ox P ± 0.003 th ± (0.0065 d5 ). (37) 

The N 4 LO estimates in (|36[> and (|37p correspond to the central value d§ = 378. The 
additional theoretical error in the coupling constant induced from the uncertainty in 
the fifth order unknown coefficient (d§ = 378 ± 378) takes the values 0.0045 (~ 1.5%) 
and 0.0065 (~ 1.9%) in the new and standard extraction procedures respectively. Com- 
paring formulas (|36p and (I37|l . one sees that within DCIPT the indicative theoretical 
errors take smaller values. In contrast to this, the experimental errors on the values of 
a s increases by the factor of 1.75 within the new procedure. It is remarkable that the 
more reliable estimate of the theoretical error presented in within CIPT (at N 4 LO) 
is close to our estimate of the error presented in formula (|37|) . 

Similarly, determining the theoretical and experimental errors on the parameter 
sj, we find stable results 

SdlNLO = 1-710 ± 0.054 cxp ± 0.002th GeV 2 
SdlN^LO = 1-709 ± 0.054 cxp ± 0.001th GeV 2 
Sdli^LO = 1-707 ± 0.054ox P ± 0.001th GeV 2 
SciIn*lo = 1-705 ± 0.054 cxp ± 0.001th GeV 2 . (38) 

It is seen from (|38[l . that the estimate for the duality point decreases very slowly with 
increasing of the order of perturbation theory. Practically, it is constant, ~ 1.71±0.05 
GeV 2 . 

Usually, it is convenient to perform evolution of the cv s results to the reference 
scale M z = 91.187 GeV. This is done by using RG equation and appropriate matching 
conditions at the heavy quark (charm and bottom) thresholds (see [48] and litera- 
ture therein). The three- loop level matching conditions in the MS scheme were de- 
rived in [49] . In this paper, we follow the work [50] ■ We perform the matching at the 
matching scale m t h = 2/% where is a scale invariant MS mass of the heavy quark 
/Lth — "Mi(Mh)- We assume for the scale invariant MS masses the estimates of the Par- 
ticle Data Group fi c = l-^ton GeV and fi h = 4.19^° GeV [ST]. In the evolution 
procedure, we use the exact numeric four-loop running coupling. In Table [3] we com- 
pare the estimates for a s (M 2 ) obtained from the two (DCIPT and CIPT) r-decay 
determinations of the coupling constant. 

As stated above, we have used the MS scheme four-loop running coupling uniformly 
in all calculations, whereas the order of approximation to the Adler function has been 
varied consecutively. To perform a more accurate test, let us now employ the same 
orders to approximate the /3 and Adler functions. The coefficient cl is accordingly 
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Table 5 Testing the stability of the results with regard to higher order perturbation theory 
corrections. Here, the approximations to the /3 and Adlcr functions are chosen consistently, at 
the same orders. 

Pcrturbative order | cl «d GeV 2 «s("i 2 )|rjciPT 

LO 0.444444 1.721 0.394 

NLO 0.336798 1.713 0.335 

N 2 LO 0.527261 1.709 0.321 

N 3 LO 0.651373 1.707 0.313 



Table 6 Comparison of the RSI and DCIPT determinations of the MS coupling constant from 
the r-decay data. Experimental errors are given only. 



Perturbative order 






NLO 

N 2 LO 

N 3 LO 


0.278 ± 0.003 
0.319 ±0.004 
0.312 ±0.004 


0.335 ±0.016 
0.321 ±0.016 
0.313 ±0.014 



recalculated. In Table we present the results of the improved test. Comparing the 
numbers in Tables [2] and [5] we see that the extracted values for the parameters, beyond 
LO, are very close (the N 2 LO and N 3 LO results practically coincide). 

Let us now employ the renormalization scheme invariant extraction method (RSI) 
of [H] to extract the numerical values of the coupling constant from the 2005 ALEPH 
V±A spectral data. We shall also include into consideration the recently calculated 
C(af ) term in the series expansion of the Adler function. The advantage of this tech- 
nique is that one starts from the physical quantity, the effective charge defined by 

a r = ^=6$ (39) 

where is the perturbative correction to the r-decay rate. The running coupling a T 
defines the internal scheme for the physical quantity. The numerical value for the QCD 
scale in the internal scheme, A T , is extracted using the equation flr(7/i^) — ^cxp- 
The MS scheme scale parameter is determined according to the relation ^jg- = 
A T exp{— 5.20232/(2/?o)}, where /3o = 9/2. Formulas for calculation of the coefficients 
of the function /3 T (the /3-function in the internal scheme) may be found in works [3] 
141] , For the experimental value of the perturbative part of the r decay rate in the 
non-strange channel, we use the updated value 

Si%\v+A = 0.2042 ± 0.0050 cxp , 

evaluated recently in [23]. For consistency, we use the same orders to approximate 
the f3 and Adler functions in the MS-scheme. In Table (0, we compare, the RSI and 
DCIPT determinations of the coupling constant order-by-order in perturbation theory. 
The relevant channels which have been used to extract the coupling are indicated by 
subscripts. It is seen from the Table, that beyond NLO the two determinations of the 
coupling constant are in good agreement. 

As is known, mathematically, the extraction of QCD parameters from experimental 
data via sum rules constitutes a so called ill posed inverse problem (analytical contin- 
uation of an approximately known function) [T7]. Small changes in the input data may 
lead to large changes in the output. In this regard, it is desirable to check the new 
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framework. To do such a test, one may extract the values of the parameters using 
the data from different r-decay experiments. As a different experimental data, let us 
employ 1998 OPAL experimental data on the non-strange isovector vector spectral 
function which is publicly available 0. The data are arranged in 100 bins with bin size 
0.032 GeV 2 , starting from s = 0.016 GeV 2 . Note that, the OPAL data correspond to 
the branching fractions available in 1998, as well as the then-current values of V n d and 
the electronic branching fraction B e . These parameters have been updated since then. 
The 2005 ALEPH analysis is more recent and based on more statistics. However, it 
was pointed out in [19] that the correlations due to unfolding have been omitted in the 
original ALEPH analysis. Thus the publicly available covariance matrices [8] should 
be corrected. Fortunately, the above mentioned obstacles have little relevance to the 
problem under investigation. In fact, our aim is to investigate the impact of the specific 
formulation of the quark-hadron duality (as given in l|18p) on the extracted value of 

Inserting into system of equations (13 2p~(|33p the empirical vector spectral function 
reconstructed from the 1998 OPAL data|^j, we solve the system numerically. We use 
the N 2 LO and N 3 LO approximations to the Adler function combined for consistency 
with the three- and four-loop order MS running couplings respectively. To determine the 
experimental uncertainties on the extracted values of the parameters we use covariance 
matrices provided by OPAL (relevant formulas were derived in the appendix to |40| ) . 
For the duality point, we obtain stable result 

Sdli^LO = (1-680 ±0.100cxp) GeV 2 

SdlN^LO = (1-679 ± 0.100e X p) GeV 2 , (40) 

the central value in (|40p decreases very slowly as the order in perturbation theory 
increases. We see that the numerical values for the duality point extracted from the 
ALEPH and OPAL data are close (cf. (|38)l and (|40)l . For the strong coupling, from the 
OPAL data, we find the values 

a s (m?)| N =LO = 0-296 ± 0.025cx P 

a s («v)lN3LO = 0.290 ± 0.023 OX p, (41) 

the central values here are somewhat smaller as compared to the corresponding values 
extracted from the ALEPH data (cf. formulas (|36p and (|4ip ~). However, the two deter- 
minations of the coupling constant are consistent within their mutual errors. It should 
be remarked that, in the case of the OPAL data, we have obtained larger experimental 
uncertainties on the numerical values of the parameters. For comparison, the original 
OPAL analysis of the same data, within CIPT, gave the value [9] 

« s (m 2 )j N 3 LO = 0.348 ± 0.009 CX p ± 0.019 th , (42) 

Comparing the numbers in formulas (|4ip and Q42p . one sees that the DCIPT deter- 
mination of a s is significantly smaller, and the two determinations are not consistent 
within their mutual errors. Performing evolution of the a s values (|41l) to the Z°-mass 
scale, we obtain 

a s (M 2 )j N 2 LO = 0.1154 ± 0.0034 cxp ±0.0005 ev 

a s (M 2 )| N3LO = 0.1145 ± 0.0033 cxp ± 0.0005 cv . (43) 

9 I would like to thank S. Menke and S. Peris for making the data available to me. 
10 We use cubic splines to interpolate the data. 
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4 Conclusion 

We have extracted numerical values for the MS scheme strong coupling constant a s 
from the r-lepton decay data. The data provided by 2005 ALEPH and 1998 OPAL 
experiments are employed. We examine in detail the dispersive approach to the r- 
decay suggested in our earlier work [40] . The errors observed in some numerical results 
of HO] have been corrected. Accordingly some of the conclusions of [3D] are changed. 

The new framework is based on the approximations to the Adler function which 
have correct analytical properties. So that the application of the FESR (O is math- 
ematically justified. Moreover, these approximations correctly reproduce the infrared 
and ultraviolet behavior of the exact Adler function. In contrast, in the standard ap- 
proaches (FOPT, CIPT or APT) some of these properties of the Adler function are 
violated. The global quark-hadron duality is used in the limited region of the energy 
squared < s < (E^ = ^/Sd ~ 1-31 GeV). In the region < s < s<j, the hadronic 
spectral function is reconstructed from the experimental data. This enabled us to re- 
duce the effects of duality violations coming from the low energy region. In fact, one 
expects in this region sizeable non-perturbative corrections to the Adler function. 

Technically, the new method is based on the system of equations (132ll -(133 p . The 
first equation follows from the concept of global quark-hadron duality employed on 
the limited interval of the energy squared, s<j < s < m 2 -. The second equation is a 
consequence of the OPE which imposes the restrictions on the ultraviolet behavior of 
the Adler function. The parameters a s and are simultaneously extracted from the 
data. We have examined numerical stability of the extracted values of the parameters 
order-by-order in perturbation theory. The new framework (DCIPT) and the standard 
(CIPT) are systematically compared. We have demonstrated that the DCIPT deter- 
minations of the strong coupling constant are more stable against perturbation theory 
corrections (see Table [3}. The central value of the coupling constant definitely became 
smaller as compared to the CIPT result (cf. Eqs. ()36[) and H37|) ~l. The changes in the 
central values are not within the quoted experimental and theoretical errors. Using the 

error estimated within DCIPT, a — <J <7c Xp + <r 2 h » 0.0151, we find that at N 3 LO the 

central values of a s (rn^) in formulas Q36p and (|37|1 differ from each other in about 2.7 
standard deviation. However, assuming the error estimated within CIPT , a w 0.0107, 
one finds even large difference, 3.8 <t[1_j. A shortcoming of the new procedure is the in- 
creased experimental error on the extracted values of a s . This is a direct consequence 
of the reduction of the duality region. 

Having included into analysis the fourth order coefficient and the geometric 
estimate d§ = 378, we have observed excellent agreement between the lattice and r- 
decay determinations of the strong coupling constant. At N 4 LO, the central value for 
a s ( see Tabled} coincides with the central value quoted in [B] (see formula @). For 
this reason, we believe that DCIPT provides better approximation as compared to 
CIPT. 

For comparison purposes, we have extracted the strong coupling constant from 
the 2005 ALEPH V+A data by using the RSI method of work gT], extending the 
result of [41] up to N 3 LO (see Table [6}. Good agreement between the RSI and DCIPT 
determinations of a s has been observed. 

The duality point is found to be surprisingly stable with respect to higher order 
QCD corrections: s d = 1.71 ± 0.05 oxp ± 0.00 th GeV 2 (see Tables [2] and . In Tabled 

11 Due to the large experimental error within DCIPT, ct(DCIPT)/ct(CIPT) w 1.4. 
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we have performed a more accurate test of stability of the numerical results, choosing 
consistently the orders of the approximations to the /3- and Adler functions. 

To examine the stability of the numerical results with respect to change in the input 
data, we have also analyzed the 1998 OPAL data for the non-strange vector spectral 
function. The extracted values for the parameters from the ALEPH and OPAL data 
are found to be consistent. 

The procedure suggested here can obviously be extended for analyzing the non- 
strange r-data from the axial- vector (A) and vector plus axial- vector (V+A) channels. 
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